The validity of the Landau-Zener model for output couphng of Bose condensates 
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We investigate the validity of the Landau-Zener model in describing the output coupling of Bose 
condensates from magnetic traps by a chirped radiofrequency field. The predictions of the model 
are compared with the numerical solutions of the Gross-Pitaevskii equation. We find a dependence 
on the chirp direction, and also quantify the role of gravitation. 
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Bose- Einstein condensation of alkali atoms in magnetic 
traps was first observed in 1995 ||l[. The first step to- 
wards coherent matter beams was taken when a con- 
trolled release of condensed atoms from a magnetic trap 
was demonstrated in MIT [gj, and recently at Garch- 
ing p]. The basic tool for output coupling has been 
spin-flipping induced by a radiofrequency (rf) magnetic 
field [^] , but other approaches have also been used [^,|j . 

The MIT group demonstrated output couplings based 
on chirping the rf field, and on resonant rf pulses M. 
In both cases the experimental results are understood in 
terms of simple models. The theory of the pulsed and 
cw output couplings has been studied extensively in the 
literature [p|-pO|. In this Brief Report we examine the 
chirped output coupler. We find the validity conditions 
for the multistate Landau-Zener (MLZ) description ||ll| ] 
and compare them with the numerical solutions of the 
Gross-Pitaevskii (GP) equation |lj]. We chart the com- 
binations of the chirp speed A = duj-cf/dt and the rf field 
amplitude Q. for which the MLZ description fails. We 
show that this failure depends on the chirp direction. 
Also, the repulsion between particles in the condensate 
and gravitation are found to contribute to the validity 
conditions. 

Magnetic trapping is achieved by spin-polarising the 
atoms. For simplicity we consider the case of F = 1, 
where F is the hyperfine quantum number. In the in- 
homogeneous magnetic field the three substates M = 
— 1, 0, 1 experience the potentials shown in Fig. |l|(a). The 
atoms on the M = ^\ state are trapped in the harmonic 
potential. In the MIT trap one had lOx = oJz = (Stt) 320 
Hz and Uy = (27r) 18 Hz ll^]. When the spin of an atom 
is flipped from this state to the Af = state, the atom 
moves away due to quantum mechanical dispersion and 
the repulsion between the condensed atoms. If the inter- 
nal state of the atom is changed into the M = 1 state, it 
feels the inverted parabolic potential as well. 

We can eliminate the rf oscillations by making the ro- 
tating wave approximation, and then shifting each M 
state in energy by an appropriate number of photons [[l4| . 
Now the field-induced resonances appear as potential 
crossings. At the rf field frequency wq the atoms at the 
center of the trap are in resonance with the field. We 
define the field detuning as A = Wj-f — wq. For A < 
none of the atoms are in resonance, and for A > atoms 
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FIG. 1. Position dependence of the energy levels of atoms 
in a magnetic trap, (a) No output coupler, (b) output coupler 
with luq > LUri, (c) output coupler with wo < Wrf- 



at locations x — ±x' are in resonance [Fig. |l|(b) and (c)]. 

In the chirped output coupler one sweeps Wj-f so that all 
atoms in the trap feel a resonant field for a brief moment. 
The idea is to make this moment long enough for achiev- 
ing a total or partial spin flip, and brief enough that the 
atoms do not have time to move due to changes in their 
internal state. Then the atoms remain stationary while 
they experience a time-dependent change of the energy 
difference between adjacent spin states. 

In the MIT experiment a linear chirp was used: A — 
Xt. For two internal states this corresponds directly to 
the Landau-Zener (LZ) model. In moderate magnetic 
flelds the Zeeman shifts are linear and thus all adjacent 
states are resonant simultaneously, so one has a genuine 
multistate problem |0,n]. In this particular case, how- 
ever, the spin dynamics of a stationary atom can be de- 
scribed analytically, by a simple multistate generalization 
of the Landau-Zener result (MLZ) |ll|. In the experi- 
ment the agreement with the MLZ prediction was good. 

For simplicity we consider one spatial direction only. 
The three-state Hamiltonian for a stationary atom lo- 
cated at X is 
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Here m is the atomic mass, lu is the trap frequency and 
hCl is the rf field coupling. For a stationary atom we 
can solve the time-dependent Schrodinger equation with 



H. The MLZ theory predicts the final populations of the 
spin states: 



P-i ==exp(-27rr), 

Po = 2exp(-27rr)[l -exp(-27rr)], 
P+i = [1 - exp(-27rr)]2, 
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where T = n'^/{2X) ||,|lT|. In the MIT experiment 
A — (27r) 500 MHz/s, and Q. changed from to about 
(27r) 11.3 kHz. The Pi's depend only on F and not on 
the trap geometry; thus in asymmetric traps the output 
coupling takes place in all directions with the same effi- 
ciency. 

In the Gross-Pitaevskii theory the single atom ampli- 
tudes ^i(a;,t) describe effectively the whole condensate 
(i.e., |^i(a;, i)p gives the density distribution of the atoms 
on spin state i) ||l3]. Their time evolution is obtained 
from the time-dependent Schrodinger equation with the 
Hamiltonian 
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The last term describes interactions between the par- 
ticles. The parameters Ci^i are proportional to atom 
numbers and the scattering lengths of the correspond- 
ing internal states i and I. In our one-dimensional study 
this parameter does not match properly with the realis- 
tic three dimensional situation. But a trap with a very 
low frequency at one direction can be regarded quasi one- 
dimensional. In this case a reasonable estimate for our 
ID C parameter would be C = Csd/^, where A is the 
cross sectional condensate area. As an order of magni- 
tude estimate it should be valid for other traps as well. 
For simplicity we take all C"s to be equal. 

We solve the GP equation numerically for the output 
coupler. In the limit of large C we can ignore the ki- 
netic energy term and obtain the Thomas-Fermi solution, 
^(a;,0) = •\/[/i — U{x)]/C, where U{x) is the trapping 
potential. The condition fj, — U{x) > defines the edge 
of the condensate. The chemical potential fi is obtained 
from the normalisation of the wave function. 

A breakdown of the MLZ model is expected if the 
atoms move during the transition process. A similar 
problem arises for diatomic molecules interacting with 
short laser pulses |jlq]. In order to quantify this break- 
down we consider the characteristic time scale 6t of the 
LZ process, St ~ fl/X |1^. The atoms need to remain 
stationary during this time. The term "stationary" can 
be defined by transforming St into a region Sx around the 
location xq of the atom. For simplicity we assume that 
xo > 0. For parabolic potentials in the P = 1 case we set 
hfl = ^muj'^[{xo + SxY — Xq] ~ muj'^XoSx, which defines 
Sx. If the atom moves a distance Ax in time St, it can 
be regarded stationary if Aa; <C Sx. 

The atomic motion can arise from quantum mechanical 
diffusion, repulsion between atoms, or from acceleration 
along the inverted parabolic potential. We consider the 



acceleration a first. With Newtonian dynamics we get 
a{xo) = —{l/m){dU/dx)\x=xo = uj'^xq. In the small re- 
gion around xq we have Ax = a{St)'^/2 = oj'^Xf)^^/{2X^). 
Thus we get the condition A^/J7 ;» mio'^x^/fi. This needs 
to be true for all a;o; the right-hand side is maximised at 
the edge of the condensate. 

For small C the edge is near the width of the ground 
state of the harmonic potential, max(a;o) — yhjrrujj, 
which gives 
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For large C we take the Thomas-Fermi approximation, 
and then max(xo) — y^^iTfrmJ^, which gives 
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where /i — huj^. The MIT trap parameters satisfy con- 
dition (6) well in all directions, and to break down con- 
dition (7) would require an unrealistically large ^. Since 
linic^o M = fiw/2, Eq. (^ is a special case of Eq. (^). 

Diffusion and repulsion can give atoms a velocity which 
initially overcomes the acceleration. For the M = state 
these processes are naturally covered by the various stud- 
ies of the ballistic expansion of condensates; see Ref. |17| ] 
and references therein. As we are only looking for con- 
straints it is sufficient to characterise the maximum speed 
V of the atoms with the energy stored in the trapped con- 
densate. We set /x = mw^/2. For small C the speed v 
reduces to the free-space momentum width of the Gaus- 
sian harmonic oscillator wave function, v ~ ^Jfiio jm. 
Now Aa; ~ vSt. On th e other hand, at xq = we 
have Sx = \j2Ti^/{mijj'^). Eventually we get for diffu- 
sion/repulsion the conditions (||) and (Q) . 

These conditions do not depend on the direction of the 
chirp. However, there exists another breakdown mecha- 
nism for the MLZ theory. Atoms that have interacted 
resonantly with the field can re-enter the resonance re- 
gion (a reunion) due to their motion. Let us consider a 
positive chirp: a resonance emerges at x = and then 
separates into two points that move towards large |a;|. 
This is demonstrated in Fig. [^ if we consider it as a se- 
quency of snapshots. Strong acceleration on the AI = 1 
state leads to the above problem. We need a condition 
on A and f2 for avoiding the reunion until the transition 
probability is negligible. 

Here the direction of the chirp is crucial. For nega- 
tive chirps the resonances emerge at large |a;l and move 
towards x — Q where they disappear. Acceleration, how- 
ever, moves the atoms in the opposite direction. Thus for 
negative chirps the reunion problem is absent. The role 
of chirp direction has been studied e.g. in the context of 
laser- molecule interactions |l5|| . 

We consider the case where the reunion happens at 
X = Xq. Energy conservation gives us roughly the local 
speed Wo of the moving atoms: wo — loxq. The slope of the 
local energy difference between adjacent states is a^ = 



muj'^Xo. The product of these quantities gives us the local 
motion-induced change in the energy levels for the atom: 
h\o = \dU/dt\ = \{dU/dx){dx/dt%=^^ = muj'^xl. The 
motion of the resonance, A, is small compared to Aq for 
realistic parameters, and we can ignore it. Basically, we 
want that the motion-induced LZ transition probability 
at xq is smaller than some fixed value F: 



r = -— < r. 

Ao 
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The time it takes for a resonance to reach xq is to = 
mu)'^XQ/{2hX). For a accelerating atoms we need to solve 
the Newtonian equation of motion: [cPx/dt^) —uP'x = 0. 
For x' ^ .J2n\t'/{muj'^), v{t') = v' we get 

v' 
x = x' cosh[w(t - t')] + — sinh[w(t - t')]. (9) 

LO 

There are several possible values for x' and v' , and the 
quantum mechanical diffusion/repulsion complicates this 
simple Newtonian picture. For large C we can assume 
that the particles which reach the reunion first come from 
the edge of the initial condensate. This fixes x' (and t'). 
We have simulated the problem numerically with the GP 
equation, and obtained the reunion time for the fastest 
atoms. As shown in Fig. 0(a) an effective trajectory v' oc 
C'^/^ locates the reunion well. The main dependence on 
C, however, arises from the location of the condensate 
edge, x' oc C^^'^. Some examples of the breakdown are 
shown in Fig. 0(b). 

One should not make detailed conclusions from such 
trajectories. Apart from the crudeness of the Newtonian 
mechanics, the fastest atoms are only a fraction of the 
condensate. In Fig. we show in the (A, J7) plane where 
the MLZ prediction for P+i fails for more than 10 % 
for the large C case. The constraint obtained by using 
the effective trajectory and condition (g) is clearly too 
demanding. Also, in a real experiment one can switch the 
field off before the reunion and thus avoid the problem. 
This is the case in the MIT experiment: after reaching 
the resonance ujq the rf field is on for about t ~ 0.5 ms. 
As TLdx — 1 and TUy = 0.06, the field is off by the time of 
reunion, as Fig. 0(a) shows. For a negative A we saw no 
large deviation from the MLZ prediction for the values 
used in Fig. 0. 

The effects of gravitation should be considered in out- 
put couplers HJl^l . In the direction of gravitation (z) the 
trapping potentials are 
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where g is the gravitational acceleration. The M = ±1 
potentials remain harmonic but have spatially shifted 




200 



0.9 
0.8 
0.7 
0.6 

0.4 
0.3 
0.2 
0.1 




















/ 


X 




MLZ 




- C/Co=10 












/ 








C/Co =100 










/ 






^ 














/ 




o 














/ 


o 






- 








/ 


O 


- 


^ 




- 






/ 


A 










- 






^1 














/ 


A 



































10 
£1/(0 



15 



20 



FIG. 2. (a) The time it takes for the spin-flipped atoms to 
reach the reunion as a function of C, and when and A = 1000;"^ . 
The circles are the GP results, and the dashed line is the 
classical estimate based on Eq. (bl). (b) The breakdown of 
the MLZ theory as a function of Q, for two values of C, and 
A = 200(.j^ Here Co 
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centers, as shown in Fig. ^ As gravitation affects only the 
external degrees of freedom the rf field resonance condi- 
tions still follow the purely magnetic potentials. In other 
words, the resonance points are located symmetrically in 
respect to the magnetic field minimum, but not in respect 
to the trap center. For a stationary atom the LZ param- 
eter F remains unaffected by gravitation, but many other 
conditions change. 

For the MIT trap one gets (Ax)g — g/oSl — 2.5 /im 
which is about the size of the trap ground state (1.7 /xm) 
but clearly smaller than the condensate (17 /im). For 
gravitation to dominate acceleration on the M = 1 state 
at the condensate edge we obtain the condition 2^uj1 < 
mg^ (for large C), which reduces to hu'^ < mg^ as C ^ 0. 
If gravitation dominates, then the basic validity condition 
for the MLZ approach becomes (Az — gfi^/A^, Sz = 
nn/{muj'^zo), zo ~ g/w^) 
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This applies also for the motion on the M = state. 
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FIG. 3. The validity of the MLZ theory. The model pre- 
dictions for P+i are compared with the GP solutions for the 
case of positive (A > 0) chirp direction. Here C = lOOCo, 
which corresponds well to the Thomas-Fermi limit. The plus 
signs (-I-) indicate when the values differ by more than 10 %, 
and circles when less. In the region below the solid line the 
reunion is not expected to produce more than 10 % devia- 
tion even if all atoms on the M = -\-l state would reach the 
reunion as soon as the first ones do it in simulations. 



Figure shows also that the shifts in the potential 
centers can be used for directed output coupling. With 
a very slow negative chirp one can leak the condensate 
from the earthside edge of the trap, as the chirped rf pulse 
becomes resonant there first. Due to the slowness the 
condensate edge will follow the resonance, and the other 
resonance point will not reach the diminishing conden- 
sate. This approach is used in the recent experiment [pj. 
It is complementary to the situation considered by us, as 
there the chirp timescale must be clearly longer than the 
motional time scales. 

In this Brief Report we have derived the validity 
conditions of the multistate Landau-Zener approach for 
chirped output couplers, with and without the presence 
of gravitation. Comparison with numerical results shows 
that these conditions allow one to achive "safe" parame- 
ter regions easily in the experiments. Our study is for one 
dimension only, but our results should apply directly to 
three dimensions, especially in the case of strongly asym- 
metric traps, where the tightest trapping direction will 
dominate the expansion of the untrapped atom cloud. 

We acknowledge the support by the Academy of Fin- 
land. 
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